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Abstract 

We derive general analytic formulae for the power spectrum and spectral index of 
the curvature perturbation produced during inflation driven by a multi-component 
inflaton field, up to the second order in the slow-roll approximation. We do not 
assume any specific properties of the potential or the metric on the scalar field 
space, except for the slow-roll condition, Einstein gravity, and the absence of any 
permanent isocurvature modes. 
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1 Introduction 



Inflation of the early universe [1] magnifies microscopic quantum fluctuations in the infla- 
ton field cj) into macroscopic classical perturbations in space-time and matter. The latter 
are supposed to be the seeds that grow to become the rich structures, such as galaxies 
or clusters of galaxies, that are observed today. Thus the power spectrum P and spectral 
index n predicted by a model of inflation can be tested from observation of the large-scale 
structures, and it is therefore important to calculate them as accurately as possible. 
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Standard calculations [2] of P and n work to lowest order in the slow-roll approximation, 
and assume that has only one dynamical degree of freedom. However, the latter as- 
sumption has no theoretical or observational justification. Previously, Stewart & Lyth [3] 
computed P and n up to the second order in the slow-roll approximation, in the single 
component infiaton case. Sasaki & Stewart [4] derived general formulae for P and n in 
the multi-component infiaton case, but only up to the first order in the slow-roll approx- 
imation. In this paper, we calculate P and n for a mult i- component infiaton up to the 
second order in the slow-roll approximation. The method of our calculation also refines 
the earher one in Ref.[4]. We use the units c — h — SnG — 1. 



2 Homogeneous background and slow-roll approximation 



Let 0" be the multiple scalar fields that slowly roll on the potential during infiation. 
We start from the action of the form: 



/ 



(1) 



where g^^,, and R are the metric and curvature scalar in space-time, and hat is the metric 
on the 0-space (which may be curved). In an exactly homogeneous universe, we have 
a^c/)" = (i = 1, 2, 3) and 

ds^ := gf.^dx^dx" = dt"^ - a'^{t)Sijdx'dx^ (2) 



(we assume that the background universe is spatially fiat). Then, from Eqs.l and 2, 
4>"' := dQ4>°'{t) obeys the equation 

0" + 3//0" + y-" = (3) 

where 0" := {D/dt)(j)°' :— ^''Vb^", Va is the covariant derivative operator associated with 
hab, and H := d/a is the Hubble parameter. Raising and lowering of indices abc- ■ ■ are 
always done by hab- Also we have from the Einstein equation 

3H'^^r<Pa + V, (4) 

H^-l^^k- (5) 

We assume that the potential 1^(0) has a sufficiently gentle slope: 

VaV''' < y^ (VabV'^^'y/'' < V , etc. (6) 
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(we list all the assumptions we need in this paper more rigorously in Appendix A). Then, 
from Eqs.3 and 4, 4>°'{t) soon approaches a slowly rolling state given by 

— ~ ~ . 7 

H 3H^ V ^ ^ 



Let us define 



These are small quantities (<^ 1) of the same order in the slow-roll approximation. They 
are slowly varying and their time derivatives (divided by H) are smaller quantities of the 
next order. For example 

H-^a = 2a{a + (5) (9) 



is a second order quantity. We use the notation ~ when the equality is valid only up to 
the lowest order in the slow-roll approximation. 



3 Perturbation 

A scalar perturbation in the space-time metric is most generally written as [5] 

ds^ = (1 + 2A)dt^ - 2a{diB)dtdx' - a^[(l + 2n)5ij + 2didjE]dx'dxK (10) 



Here TZ is interpreted as the intrinsic-curvature perturbation in the constant time hyper- 
surface. Let (50" be the perturbation in the scalar fields around 0"(t). In appendix B, we 
derived from Eqs.l and 10 the equation of motion of (50" on fiat hypersurfaces: 



dt-^ at \a a"^ dt 



a lb 



H 



0>' 



where R'^bcd is the Riemannian curvature tensor in the 0-space. We work in fc-space 
throughout and simply use (50" as the Fourier transform of the perturbation. The conformal 
time T) is defined by 

f dt 1 f ada 
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Since a is slowly varying (see Eq.9), we take a out of the integral and obtain 



r}= -{l + a)/{aH). (13) 

Defining u°' :— a&jf and working to the first order in a, /3, etc., Eq.ll is rewritten as 

+ k\'' = ^(2m« + 3ey) (14) 



where we regard 

^ab '■— Oiriab + [llacllbd — g-^acbd I -Jp [^^) 



as a first order quantity. In order to solve the differential equation 14, we introduce 
the orthonormal basis e'^ {A runs over the number of scalar field components) parallel- 
transported along the unperturbed trajectory 0"(t): 

^.-0 (16) 



SO that the symmetric tensor is diagonalized along 0"(i) as 

e^'^Y.'Ael^e",). (17) 



This diagonalization is justified as follows. At some point on 0"(i), can be diagonalized 
exactly as in Eq.17 with the eigen- vectors e". As one moves along 0"(i) with parallel- 
transported, will not remain the eigen- vectors and off-diagonal components may appear 
in Eq.17. However, since we are assuming (see Eq.A2) that Cab is covariantly changing 
slowly along 0°'(t), the off-diagonal components are second order quantities. Therefore, 
Eq.17 is valid up to the lowest order in the slow-roll approximation. In short, we treated 
as a constant in Eq.17, just as we treated a as a constant in Eq.l2. Prom Eqs.16 and 
17, the 74-component of Eq.l4 is written as 

^ + ^ ^{2 + 3e,)u, (18) 

where := «ae"- First let us consider microscopic fluctuations, the physical wavelength 
of which is well-inside the horizon {—krj 00). When —kr) 00, the right hand side 



4 



(RHS) of Eq.l8 is negligible compared with the A;^ term, and thus the Ua behave like 
(real) massless Klein-Gordon fields: 



(19) 



where a\ and are the creation and annihilation operators of an A-particle: 

[a^(fc), ai{k')] = d^sS^k - k'), ajO) = 0. (20) 



As —kf] approaches unity, the RHS of Eq.l8 becomes comparable to the A;^ term, and the 
solution is written in terms of the Hankel functions as 



(fc) = i-kriy/' \CA{k)Hl^l{-kv) + Ci{-k)Hi^^l{-kri) 



(21) 



where 



(22) 



Using the asymptotic behavior of the Hankel functions at infinity, the integral constants 
Ca are determined from Eq.l9 as 



c^(fe) 



exp 



{2ua + 1) 



aA{k) 



2k 



(23) 



Next we go well-outside the horizon, i.e., —kr] [but — \n{—kr]) is not too large]. Using 
H^^''^\x) ±r{i/){2/xY /{m) as X ^ 0, and expanding up to the first order in e^, one 
finds 



^ ^ V2k\k7]) r(3/2) \ kr]J ^ ^ ^ 



' (^]{l + [c-ln{-kr))]eA}bA{k) (25) 



V2k \kri 
where 

bA{k) := e^"^^/^a^(fc) - e-^^^^^/^a^-fc) , (26) 



and c:=2 — ln2 — 7 = 0.7296 • • • with the Euler number 7. It is clear from Eq.25 that 
the perturbations become completely classical as we go outside the horizon, because 
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and a\ appear only in the combination of 6^ and hence [li^jTiyi] = follows. Going back 
to general coordinates, we obtain 



c + In 



/aH 



(27) 



with 6" := bAe'^. For later use, we note that 

{ba{k)bl{k')) = habS'ik - k') 



(28) 



where (• • •) reads the vacuum expectation value. 



4 Power spectrum and spectral index of TZc 



Sasaki & Stewart [4] showed that the curvature perturbation on a comoving hypersurface 
TZc (see Eq.lO) during the radiation-dominated phase (after complete reheating) is related 

to (5^" as 

7^e = iya<5i/'^ (29) 



The RHS is to be evaluated at some time (say ti) during inflation soon after the scale of 
the perturbation goes well-outside the horizon (but does not depend on the exact value 
of ti, as shown below). It is also assumed in Eq.29 that the space-time is foliated on a flat 
hypersurface at ti, in accord with Eq.ll. Here 

N{(P) := I Hdt (30) 



is the number of e-folds in the background universe, and t2 is the time corresponding to 
some fixed energy density during the radiation-dominated phase. In general N can depend 
on both (j)°'{ti) and 0"(ti). However, as we are assuming that slow-roll has been achieved, 
the 0-dependence should be eliminated using the slow-roll trajectory which is given in 
Eq.C12 up to second order. The power spectrum P{k) of TZc is defined by 

{n,(k)nl{k')) = 27r^k-^P{k)5^{k - k') . (31) 



Prom Eqs.27, 28 and 29, one finds 



eabM' 



ab 



(32) 
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where Mab '■— NaNb/N''^Nc- Thus the spectral index is 

dlnP 



n := 1 + 



dink 



1 - 2e„,M'^^ 



(33) 



which is the identical result with Ref. [4]. In order to calculate n up to the second order, 
we rewrite Eq.32 to show more explicitly that P{k) does not depend on time. We expand 
the pref actor of Eq.32 as 



1 + 



dlniN'^NcH^ 



o?lna 



In 



aH=k 



(34) 



Substituting Eq.34 into 32, and using Eq.Cll, we see that the \n{aH/k) terms cancel and 
obtain 



H 



P(A;) = iV'-^iV, ( — ) (l-2a + 2ce„;,M»'') 



(35) 



aH=k 



In this expression for P{k), the /c-dependence of the LHS is such that the RHS (which 
is a function only of time) is evaluated at the horizon-crossing time aH = k. Thus P{k) 
does not depend on time, as noted above. To avoid any confusion, let us define Q{a) to 
be the RHS of Eq.35 so that 



P{k) = Q{a)\aH=k 



(36) 



Then, using Eqs.9, CIO and C14, n is calculated up to second order as 



n = 1 + — — = 1 + (1 + a)- 



dliaaH 



din a 



1 - 2q; + 2XabM''^ - 2(3 - 2c)a' - 4(1 - c)ap + -aXabM"" + 4c(AabM' 
--(6c - l)M„,AAc - 3]^^^ -J^r 3^'' + '^^dt^^' 



ab\2 



(37) 



(38) 



where 



= ahab - eab ■ (40) 

[Note that l/{N'°'Na) is a first order quantity.] The /c-dependence of n is understood in 
the same way as in Eq.35. Rewriting Eq.38 in terms of V, we find 
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where 



U :=\nV, 



(42) 



(43) 



{UaUb and U-ab are first order quantities.) 



5 Summary 

We have derived general analytic formulae for the power spectrum P (Eq.35) and spectral 
index n (Eq.38 or 41) of the curvature perturbation TZc produced during inflation driven 
by a multi-component infiaton field, up to the second order in the slow-roll approximation. 
Once one specifies a model of inflation and calculates the number of e-folds N{(f)) (Eq.30) 
in the background universe, then the substitution of them into our general formulae 
immediately yields the power spectrum and its index with accuracy. We anticipate that, 
to lowest order, Na can be calculated from the inflationary phase only. However, to the 
next order, as considered in this paper, the contribution to Na from the reheating and 
radiation-dominated phases should be signiflcant. 

The magnitude of the first order terms is of order A^^^^id in many infiation models {Ncnd 
is the number of e-folds from the horizon-crossing time to the end of infiation), and if 
thermal infiation [6] occurs after ordinary infiation, we have A'end ~ 30-40. In this case, 
the correction terms in Eq.35 should be observable, while those in Eqs.38 and 41 may be 
marginally observable. Thus, our formulae may be useful in testing models of infiation 
when P and n are observed accurately by presently-planned experiments such as the 
Microwave Anisotropy Probe. At the same time, the theories and models of infiation need 
to progress to make more precise predictions. 

We would like to thank N. Sugiyama for helpful discussions, and D. H. Lyth and Y. Suto 
for useful comments. TTN and EDS are supported by JSPS Fellowships at UTAP and 
RESCEU, respectively. This work is supported in part by Monbusho Grant-in- Aid for 



JSPS Fellows No. 95209. 
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Appendices 



A Slow-roll conditions 



Here we summarize all the conditions assumed in our calculation, without derivation (One 
can derive them by differentiating Eq.3 several times and using the resultant equations 
recursively). We assumed that i) the potential has sufficiently gentle slope, ii) slow roll 
has been achieved, and iii) the curvature of 0-space is not too large and is slowly varying 
[ i) is a necessary condition for ii)]. To write the conditions quantitatively, let us define 
the "norm" 



(Al) 



of a vector X"' and a symmetric second-rank tensor Yab = Yba, and introduce a small 
quantity £ <^ 1. Then, in the calculation of P, we assumed [conditions i) and iii)] 



\UaUb\\ < e, \\Uab\\ < e, WRachdU'^U^'^W < £, 
mabcU'^W < e\ macbdU^'U^XU^'^W < 



(A2) 



so that Eq.35 is valid up to order e. Also we need the extra but rather weak condition 

1 : 



< e 



5/2 



(A3) 



because of ii). The second line in Eq.A2 is necessary to make sure that the first order 
quantities are slowly varying. Similarly, in order that Eqs.38 and 41 are valid up to order 
the calculation of n assumes 

\\U.abcdU''U'''\\ < £^ \\U.abcU''''U4 < £^ \\[{RacbdU''U''').,eU'%fU'f\\ < 8^ (A4) 



in addition to Eq.A2, and 



1 



< s 



7/2 



(A5) 



instead of Eq.A3. 
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B Derivation of Eq.ll 



Prom Eq.l, the stress tensor is given by 



(Bl) 



and the Euler-Lagrange equation is 

D 



habg^''d,(f>'' + Va^O 



(B2) 



which yields Eq.3. Choosing the gauge TZ — B — in Eq.lO, the relevant components of 
the metric perturbation are (we work in fc-space) 



Sgoo = -5^'° = 2^ , Sgoi = dg""' = , 6 In^ = A-k'E. 



Oi 



(B3) 



Wc perturb covariantly with respect to 0" and define the covariant perturbation op- 
erator 6 := 3(p°''Va; for example, 6hab = 0. Since Sep"" is Lie-transported along 4>°'{t), i.e., 
[5/), 0]'' = 0, it follows that 



(B4) 



Using Eq.D.7 in Ref.[5] to calculate the perturbation in the Einstein tensor, 6Gq = 6Tq 
and SG^i = ST^i give 



2HA = 



(B5) 
(B6) 



respectively. Using Eq.3, one finds 



(B7) 



Perturbing Eq.B2 covariantly gives 
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+ {A-k^E)^a+(^^^ &i>a + &i>'V.,ab = 0. (B8) 

Here the second term is calculated as 

- ^j) = {&p'V,^'% - ^''V,5cl>'Vb)k (B9) 



dt dt 



= 25/)''0^0„;[c6] + (50^0'=;60a;c - «';c0a;6 (BIO) 
^Racbd4>'4>%'- (Bll) 



on account of Eq.B4. Therefore 



dt 



= + A;'^)0„ + 2A{(t)a + Sii-^a) . (B12) 

Prom Eqs.B6 and B7, it is easy to show that the RHS of Eq.B12 is equivalent to that of 
Eq.ll. 



C Some useful formulae 

To evaluate Eqs.34 and 37, we calculate 

din a m N'<^N^H ^ ^ 

where Na is calculated as 

A^a = 0'VfeV„iV (C2) 

- Vai^'VbN) - {Va^'){VbN) (C3) 

(C4) 

Taking the gradient of Eqs.3 and 4, one obtains 

Ha^-^{Va + ^%a) (C5) 
ka = -^[^kH,a + Ro.U^"^" + V,ab + ¥;aKc + {k;a)] ■ (C6) 
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The second term in Eq.C6 arises when one commutes the covariant derivatives of ^'^ViVc^b- 
From these equations, (pa-h — 4'b:a holds to the lowest order and thus the second term in 
Eq.C5 becomes 4>a- Using Eq.3 again, one finds 

Ha^—^a. (C7) 



We use = when the equation is valid up to the next order. Defining the first order quantity 

lab '■— [ -^'I'ac'l'bd - g-Kocfcd ) -Jp -^Jp , 



Eq.C6 can be written iteratively as 



k;a = H^ab ' \Hlllbc " ■ (C9) 

Prom Eqs.Cl, C4, C7, C9 and 39, one obtains 

~-2e„,M'^^ (Cll) 
Also it is easy to show that 

^ ^ -([/■« + \w^Ub) ~ -U'" , (C12) 
H 3 

H-^Na^\abN^\ (C13) 

H-^Ma, ~ 2{X\aM,,, - M^.M^X'^) . (C14) 
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